Dynamic imbalance detection apparatus and method using linear time-varying angular velocity model

ABSTRACT

Provided is a dynamic balancing apparatus using a linear time-varying angular velocity model, which includes a rotational shaft on which a rotational body having a rotating unbalanced mass is installed; a linear time-varying angular velocity generator which allow the rotational shaft to be rotated at a linear time-varying angular velocity; a support for supporting both sides of the rotational shaft; and a transducer for measuring applied force or vibration transferred to the support by revolution of the rotational shaft having a linear time-varying angular velocity.

TECHNICAL FIELD

The present invention relates to a dynamic balancing apparatus and effective balancing methods using a linear time-varying angular velocity model, which precisely evaluate the angular position and magnitude of an unbalanced mass of a rotational body which is rotated at the linear time-varying angular velocity.

BACKGROUND ART

All components of rotational machines are primarily manufactured along the central rotation axis of dedicated machines thereof. A dynamic balancing apparatus has been used to measure the angular position and magnitude of an unbalanced mass which exits with respect to the actual running axis after the components of the rotational machine are assembled.

A dynamic balancing method is a theoretical method which is used in the dynamic balancing apparatus to measure the angular position and magnitude of the unbalanced mass.

In order to implement the dynamic balancing, most commercially used dynamic balancing apparatuses measure only a centrifugal force component generated by the unbalanced mass when a balancing test piece is rotated at the constant speed, or indirectly measure vibration of a support generated by the centrifugal force, thereby estimating the angular position and magnitude of a rotating unbalanced mass.

The existing balancing apparatus needs an accelerating time until reaching to a certain rotating speed from the stopped state, a real balancing time which is performed at the certain rotating speed and a decelerating time which is necessary to stop the rotating state from the certain rotating speed.

Due to the accelerating time, the real balancing time and the decelerating time, there is a problem that the existing balancing apparatus is limited in productivity, i.e., the number of balancing operations per unit time.

DISCLOSURE OF THE INVENTION

It is an object of the present invention to provide a dynamic balancing apparatus and effective methods using a linear time-varying angular velocity model, which can perform a balancing operation over all time regions of the accelerating time, the constant speed balancing time and the decelerating time which are accompanied with the balancing operation of a rotational body.

It is another object of the present invention to provide a dynamic balancing apparatus and effective methods using a linear time-varying angular velocity model, which remarkably improve the number of balancing operations per unit time.

To solve the above-mentioned objects, the present invention provides a dynamic balancing apparatus using a linear time-varying angular velocity model, including a rotational shaft 10 on which a rotational body 1 having an unbalanced mass is installed; a linear time-varying angular velocity generator 20 which is installed on the rotational shaft and allows the rotational shaft 10 to be rotated at a linear time-varying angular velocity; a support 30 for supporting both sides of the rotational shaft 10; and transducers 40 which is installed on the support and measures applied force or vibration transferred to the support 30 by revolution of the rotational shaft 10 having a linear time-varying angular velocity.

Preferably, the dynamic balancing apparatus further includes a reference pulse generator 50 which generates a reference position pulse of the rotational shaft 10 so as to measure a position of the rotational body 1 with respect to the linear time-varying angular velocity of the rotational shaft 10.

Also, the present invention provides a dynamic balancing method using a linear time-varying angular velocity model, including: rotating a rotational shaft 10, on which a rotational body 1 having a rotating unbalanced mass is installed, at a linear time-varying angular velocity; measuring vibration generated by tangential force, centrifugal force, or the two forces generated by the rotating unbalanced mass of the rotational body 1 from a support 30 which supports the rotational shaft 10 rotated at the linear time-varying angular velocity; estimating frequency components of a measurement signal corresponding to the measured tangential and centrifugal forces; and estimating the angular position and the magnitude of the unbalanced mass of the rotational body 1 using components of the measured normal force and centrifugal forces and characteristics of the Fourier coefficients of the two forces.

Preferably, the angular position and magnitude of the unbalanced mass of the rotational body 1 is determined by the following equations:

$\begin{matrix} {{F_{X}\left( {\theta(t)} \right)} = {{{- F_{t}} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} + {F_{n} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}}}} \\ {= {{{mr} \cdot b \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} + {{mr} \cdot \left\{ {a^{2} + {2\;{b \cdot {\theta(t)}}}} \right\} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}}}} \\ {= {{mr} \cdot b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\cos\left( {\theta(t)} \right\}}} +} \right.}} \right.}} \\ {{mr} \cdot a^{2} \cdot \left\{ {{{\cos\left( \theta_{0} \right)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\sin(\theta)}}} \right\}}} +} \right.} \\ {{{mr} \cdot 2}\;{b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\sin\left( {\theta(t)} \right\}}}} \right.}} \right.}} \end{matrix}$ $\begin{matrix} {{F_{Y}\left( {\theta(t)} \right)} = {{F_{t} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}} + {F_{n} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} -}} \\ {= {{{mr} \cdot b \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}} + {{mr} \cdot \left\{ {a^{2} + {2\;{b \cdot {\theta(t)}}}} \right\} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} -}} \\ {= {{mr} \cdot b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\sin\left( {\theta(t)} \right\}}} +} \right.}} \right.}} \\ {{mr} \cdot a^{2} \cdot \left\{ {{{\cos\left( \theta_{0} \right)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\cos(\theta)}}} \right\}}} +} \right.} \\ {{{mr} \cdot 2}\;{b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\cos\left( {\theta(t)} \right\}}}} \right.}} \right.}} \end{matrix}$ wherein F_(x)(t) is the horizontal component of each of the tangential force F_(t) and the centrifugal force F_(n) generated at the support 30 by the rotating unbalanced mass of the rotational body 1; F_(y)(t) is the vertical component of each of the tangential force F_(t) and the centrifugal force F_(n) generated at the support 30 by the rotating unbalanced mass of the rotational body 1; m[kg] is the rotating unbalanced mass of the rotational body; r[m] is the distance of the rotating unbalanced mass of the rotational body, which is spaced apart from the center of rotational shaft 10; a[radians/s] is an initial rotational velocity; b[radians/s²] is a constant angular acceleration at the initial rotational velocity; θ(t) is the angular displacement of the rotational shaft; and θ₀[radians] is an angle from a rotation reference point of the rotational shaft.

Preferably, an unbalanced amount (m·r) due to the rotating unbalanced mass of the rotational body 1 is determined by the following equation using the tangential force and the centrifugal force, which are measured from the support and acting in the horizontal and vertical directions,

${\begin{bmatrix} {C_{F_{X}}(k)} \\ {S_{F_{X}}(k)} \\ {C_{F_{Y}}(k)} \\ {S_{F_{Y}}(k)} \end{bmatrix} = {\begin{bmatrix} 0 & {{- 2}\;{b \cdot {S_{C,\theta}(k)}}} \\ {2\;{b \cdot {C_{S,\theta}(k)}}} & 0 \\ {2\;{b \cdot {S_{C,\theta}(k)}}} & 0 \\ 0 & {2\;{b \cdot {C_{S,\theta}(k)}}} \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}},{{{for}\mspace{14mu} k} = 1},2,{3\ldots}$ wherein C_(Fx) is cosine components corresponding to Fourier coefficients of the horizontal components of the tangential and centrifugal forces acting on the support 30 by the rotating unbalanced mass of the rotational body 1, and S_(Fx) is sine components corresponding to Fourier coefficients of the horizontal components of the tangential and centrifugal forces acting on the support 30 by the rotating unbalanced mass of the rotational body 1, and these components are determined by the following equation:

${{C_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and k is an integer indicating the order of Fourier coefficients, C_(Fy) is the cosine components corresponding to Fourier coefficients of the vertical components of the tangential and centrifugal forces generated at the support 30 by the rotating unbalanced mass of the rotational body 1, and S_(Fy) is the sine components corresponding to Fourier coefficients of the vertical components of the tangential and centrifugal forces generated at the support 30 by the rotating unbalanced mass of the rotational body 1, and these components are determined by the following equation:

${{C_{F_{Y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{Y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and U_(cos) is the cosine component of the unbalanced amount (U_(cos)=mr·cos(θ₀)) and U_(sin) is the sine component of the unbalanced amount (U_(sin)=mr·sin(θ₀)),

${{C_{c,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${{S_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${C_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ {{\begin{matrix} {{1/4},\mspace{14mu}{{{for}\mspace{14mu} k} = 1}} \\ {{- \frac{k}{k^{2} - 1}},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,\ldots} \end{matrix}{S_{c,\theta}(k)}} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ \begin{matrix} {{{3\;{k/4}} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 0},1} \\ {\frac{1}{k^{2} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,{\ldots.}} \end{matrix} \right.}} \right.}$

Preferably, a sampling method of measurement signals is chosen to measure applied force signals on the support using a reference pulse generator or a rotary encoder which generates angular position pulses having a constant width (or equi-angle revolution) whenever the rotational shaft is rotated.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A and 1B are views showing components of applied forces exerted to the support and acceleration components of an unbalanced mass.

FIGS. 2 to 4 are views showing signal sampling methods using a rotary encoder or a reference pulse generator.

FIGS. 5A, 5B and 5C are views showing vibration models of the support.

FIG. 6 is a schematic view showing the configuration of a dynamic balancing apparatus using a linear time-varying angular velocity model in accordance with an embodiment of the present invention.

DETAILED DESCRIPTION OF MAIN ELEMENTS

 1: rotational body 10: rotational shaft 20: linear time-varying angular velocity generator 30: support 40: transducer 50: reference pulse generator or encoder

BEST MODES FOR CARRYING OUT THE INVENTION

Practical and presently preferred embodiments of the present invention are illustrative as shown in the following embodiment.

However, it will be appreciated that those skilled in the art, on consideration of this disclosure, may make modifications and improvements within the spirit and scope of the present invention.

FIGS. 1A and 1B are views showing components of applied forces exerted to the support and acceleration components of an unbalanced mass, and FIGS. 2 to 4 are views showing signal sampling methods using a rotary encoder or a reference pulse generator. FIGS. 5A, 5B and 5C are views showing vibration models of the support and FIG. 6 is a schematic view showing the configuration of a dynamic balancing apparatus using a linear time-varying angular velocity model.

If a rotational body is rotated at an initial rotational velocity a [radians/s] and a constant angular velocity b [radians/s²] while passing an angular reference point at a time, an angular displacement θ(t) of a rotational shaft, an angular velocity ω(t) and an angular acceleration α(t) can be expressed by the following equations. θ(t)=a·t+0.5·b·t ², ω(t)=dθ(t)/dt=a+b·t, α(t)=d ²θ(t)/dt ² =b  Eq. 1

As expressed in the equation 1, when the angular acceleration b is constant, the angular velocity ω has a linear relationship with a time. A method of realizing the balancing of the rotational body under the rotational motion, i.e., the linear time-varying angular velocity, is as follows.

As shown in FIG. 1A, mass m[kg] spaced by a distance r[m] from a rotational shaft and positioned at the angle θ₀[radians] from an angular reference point (a dot line illustrated in FIG. 1A) is moved at an acceleration a_(t)(t) in a tangential direction and an acceleration a_(n)(t) in a normal direction. As shown in FIG. 1B, these two acceleration components act as repulsive force such as tangential force F_(t)(t) and centrifugal force F_(n)(t) with respect to a support, respectively. F _(t)(t)=−m·a _(t)(t)=−mr·α(t)=−mr·b, F _(n)(t)=−m·a _(n)(t)=mr·ω ²(t)=mr·{a ²+2b·θ(t)}  Eq. 2 Explaining the numerical marks −m·a_(n)(t)mr·ω²(t) of F_(n)(t)=−m·a_(n)(t)=mr·ω²(t) with reference to the vector direction defined in FIG. 1A), since the centrifugal force F_(n)(t) and the acceleration a_(n)(t) direction act in the opposite direction, it is expressed in a Scalar expression using an negative (−) mark. In the coordinate system shown in FIG. 1A), the vector direction of the angular displacement, the angular velocity, and the angular acceleration describing the rotational movement of the mass m[kg] on the XY place should have only a component of the Z-axis direction. the vector direction of the angular displacement the angular velocity and the angular acceleration describing the rotational movement of the mass m[kg] on the XY plane should have only a component of the Z-axis direction. When the angular displacement, the angular velocity, and the angular acceleration have positive (+) values in a counterclockwise direction, they are directed to a positive (+) Z direction, and when the angular displacement, the angular velocity, and the angular acceleration have negative (−) values in a clockwise direction, they are directed to a negative (−) Z direction. With the angular kinematics, it should be appreciated that the angular displacement, the angular velocity, and the angular acceleration have directivity, and the rotation body has a perpendicular direction to a rotation surface around which it moves. In two dimensions, the Scalar expression a_(n)r·ω² is sufficient, whereas, in three dimensions, a vector expression is further required. The angular velocity, which is the vector, is in parallel with the rotational shaft. As described above, the angular velocity in FIG. 1A) is directed to a direction perpendicular to a ground. A direction to be rotated is determined by the Right Hand Rule. The direction of the angular velocity is the direction of the thumb when right hand fingers hold the rotational shaft along the arc of the circumference of the rotational shaft in a direction that the mass m[kg] is rotated and the thumb is raised to be perpendicular to other fingers.

When a velocity of the rotational shaft is increased at a constant angular acceleration, it can be understood that the tangential force is proportional to a magnitude of the angular acceleration and the centrifugal force is proportional to the sum of the square of initial angular velocity plus the angular acceleration multiplied by the angular displacement. The tangential force and the centrifugal force are respectively divided into a horizontal component F_(X)(t) and a vertical component F_(Y)(t) and then expressed as follows.

$\begin{matrix} \begin{matrix} {{F_{X}\left( {\theta(t)} \right)} = {{{- F_{t}} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} + {F_{n} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}}}} \\ {= {{{mr} \cdot b \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} +}} \\ {{mr} \cdot \left\{ {a^{2} + {2\;{b \cdot {\theta(t)}}}} \right\} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}} \\ {= {{mr} \cdot b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\cos\left( {\theta(t)} \right\}}} +} \right.}} \right.}} \\ {{mr} \cdot a^{2} \cdot \left\{ {{{\cos\left( \theta_{0} \right)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\sin(\theta)}}} \right\}}} +} \right.} \\ {{{mr} \cdot 2}\;{b \cdot \left\{ {{{\cos\left( \theta_{0} \right)} \cdot \theta}{(t) \cdot}} \right.}} \\ {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\sin\left( {\theta(t)} \right\}}}} \right.} \end{matrix} & {{Eq}.\mspace{14mu} 3} \\ \begin{matrix} {{F_{Y}\left( {\theta(t)} \right)} = {{F_{t} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}} + {F_{n} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}}}} \\ {= {{{- {mr}} \cdot b \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}} +}} \\ {{mr} \cdot \left\{ {a^{2} + {2\;{b \cdot {\theta(t)}}}} \right\} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} \\ {= {{- {mr}} \cdot b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\sin\left( {\theta(t)} \right\}}} +} \right.}} \right.}} \\ {{mr} \cdot a^{2} \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\cos\left( {\theta(t)} \right\}}} +} \right.}} \right.} \\ {{{mr} \cdot 2}\;{b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\sin\left( {{\theta(t)} +} \right.}} \right.}} \\ {{\sin\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\cos\left( {\theta(t)} \right\}}} \end{matrix} & {{Eq}.\mspace{14mu} 4} \end{matrix}$

As described above, it can be understood that these horizontal and vertical forces are changed according to an unbalanced amount (mr: an amount defined as a multiplication of a mass and a radius) and an angular position θ₀. This fact means a method of estimating the unbalanced amount and the angular position from components of the forces applied to the support 30 using an angular velocity motion which is linearly varied with respect to time, i.e., a new principle of balancing method. The balancing method of the present invention uses the components of tangential and centrifugal forces all together, and also uses characteristics of a periodic function of the two forces, i.e., characteristics of frequency components in estimating the unbalanced amount and the angular position. In other words, the applied forces having the horizontal and vertical components is comprised of a combination of the periodic functions cos(θ(t)), sin(θ(t)), θ(t)×cos(θ(t)) and θ(t)×sin(θ(t)) of the rotational angle. The periodic functions are expressed by Fourier coefficients showing the mathematical relationships of the periodic functions as follows.

$\begin{matrix} {{{{C_{c}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{1}{2}}},{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}}{{{S_{s}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{1}{2}}},{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}}{{{C_{c,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}}{{{S_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}}{{C_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ {{\begin{matrix} {{{- 1}/4},\mspace{14mu}{{{for}\mspace{14mu} k} = 1}} \\ {{- \frac{k}{k^{2} - 1}},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,\ldots} \end{matrix}{S_{c,\theta}(k)}} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ \begin{matrix} {{{3\;{k/4}} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 0},1} \\ {\frac{1}{k^{2} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,\ldots} \end{matrix} \right.}} \right.}}} & {{Eq}.\mspace{14mu} 5} \end{matrix}$

Therefore, it can be understood that resultant values obtained by integrating the components of the horizontal and vertical forces measured upon every revolution of the rotational shaft can be expressed by a combination of the above-described Fourier coefficients.

First of all, a cosine component C_(Fx) and a sine component S_(Fx) corresponding to the Fourier coefficients of the horizontal force are as follows.

$\begin{matrix} {{{C_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}}} & {{Eq}.\mspace{14mu} 6} \end{matrix}$

In the equation, k is an integer indicating an order of the Fourier coefficients.

A cosine component C_(Fy) and a sine component S_(Fy) corresponding to the Fourier coefficients of the vertical force are as follows:

$\begin{matrix} {{{C_{F_{Y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{Y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}}} & {{Eq}.\mspace{14mu} 7} \end{matrix}$

In case of k=0, that is, a mean value of the horizontal and vertical force components, while the rotational shaft is rotated once, is as follow.

$\begin{matrix} {\begin{matrix} {{C_{F_{X}}(0)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{F_{X}\left( {\theta(t)} \right)}\ {\mathbb{d}\theta}}}}} \\ {= {{{- {mr}} \cdot 2}\;{b \cdot {S_{C,\theta}(0)} \cdot {\sin\left( \theta_{0} \right)}}}} \\ {= {{{mr} \cdot 2}\;{b \cdot {\sin\left( \theta_{0} \right)}}}} \end{matrix}\begin{matrix} {{C_{F_{Y}}(0)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{F_{Y}\left( {\theta(t)} \right)}\ {\mathbb{d}\theta}}}}} \\ {= {{{mr} \cdot 2}\;{b \cdot {S_{C,\theta}(0)} \cdot {\cos\left( \theta_{0} \right)}}}} \\ {= {{{- {mr}} \cdot 2}\;{b \cdot {\cos\left( \theta_{0} \right)}}}} \end{matrix}{{S_{F_{X}}(0)} = {{S_{F_{Y}}(0)} = 0}}} & {{Eq}.\mspace{14mu} 8} \end{matrix}$

Therefore, it is possible to estimate the angular position and magnitude of the unbalanced mass from the mean value of the horizontal and vertical forces at each revolution of the rotational shaft. A cosine component U_(cos)=mr·cos(θ₀) and a sine component U_(sin)=mr·sin(θ₀) of the unbalanced amount mr are obtained by the following equation using an acceleration value b at each revolution. From the cosine and sine components of the unbalanced amount, the magnitude and angular position of the unbalanced mass can be finally expressed as follows.

$\begin{matrix} {{{{{mr} \cdot {\cos\left( \theta_{0} \right)}} = {{{- {C_{F_{Y}}(0)}}/2}\; b}},{{{mr} \cdot {\sin\left( \theta_{0} \right)}} = {{{C_{F_{X}}(0)}/2}\; b}}}{{{mr} = {\frac{1}{2\; b}\sqrt{{C_{F_{X}}(0)}^{2} + {C_{F_{Y}}(0)}}}},{\theta_{0} = {\tan^{- 1}\left( \frac{C_{F_{X}}(0)}{- {C_{F_{Y}}(0)}} \right)}}}} & {{Eq}.\mspace{14mu} 9} \end{matrix}$

In case of the same order as the number of revolutions, i.e., k=1, values corresponding to the cosine and sine components of the horizontal and vertical forces are expressed by the following equations.

$\begin{matrix} {\begin{matrix} {{C_{F_{X}}(1)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos(\theta)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{mr} \cdot b \cdot {C_{C}(1)} \cdot {\sin\left( \theta_{0} \right)}} +}} \\ {{{mr} \cdot a^{2} \cdot {C_{C}(1)} \cdot {\cos\left( \theta_{0} \right)}} +} \\ {{{{mr} \cdot 2}\;{b \cdot {C_{C,\theta}(1)} \cdot {\cos\left( \theta_{0} \right)}}} -} \\ {{{mr} \cdot 2}\;{b \cdot {S_{C,\theta}(1)} \cdot {\sin\left( \theta_{0} \right)}}} \end{matrix}\begin{matrix} {{S_{F_{X}}(1)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin(\theta)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{mr} \cdot b \cdot {S_{S}(1)} \cdot {\cos\left( \theta_{0} \right)}} - {{mr} \cdot a^{2} \cdot {S_{S}(1)} \cdot {\sin\left( \theta_{0} \right)}} +}} \\ {{{{mr} \cdot 2}\;{b \cdot {C_{S,\theta}(1)} \cdot {\cos\left( \theta_{0} \right)}}} -} \\ {{{mr} \cdot 2}\;{b \cdot {S_{S,\theta}(1)} \cdot {\sin\left( \theta_{0} \right)}}} \end{matrix}\begin{matrix} {{C_{F_{Y}}(1)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos(\theta)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{- {mr}} \cdot b \cdot {C_{C}(1)} \cdot {\cos\left( \theta_{0} \right)}} +}} \\ {{{mr} \cdot a^{2} \cdot {C_{C}(1)} \cdot {\sin\left( \theta_{0} \right)}} +} \\ {{{{mr} \cdot 2}\;{b \cdot {S_{C,\theta}(1)} \cdot {\cos\left( \theta_{0} \right)}}} +} \\ {{{mr} \cdot 2}\;{b \cdot {C_{C,\theta}(1)} \cdot {\sin\left( \theta_{0} \right)}}} \end{matrix}\begin{matrix} {{S_{F_{Y}}(1)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin(\theta)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{mr} \cdot b \cdot {S_{S}(1)} \cdot {\sin\left( \theta_{0} \right)}} + {{mr} \cdot a^{2} \cdot {S_{S}(1)} \cdot {\cos\left( \theta_{0} \right)}} +}} \\ {{{{mr} \cdot 2}\;{b \cdot {S_{S,\theta}(1)} \cdot {\cos\left( \theta_{0} \right)}}} +} \\ {{{mr} \cdot 2}\;{b \cdot {C_{S,\theta}(1)} \cdot {\sin\left( \theta_{0} \right)}}} \end{matrix}} & {{Eq}.\mspace{14mu} 10} \end{matrix}$

In the above equations, the Fourier coefficients which were presented in the equation 5 are as follows.

$\begin{matrix} {{{C_{C,\theta}(1)} = {{S_{S,\;\theta}(1)} = \frac{\pi}{2}}},{{S_{C,\theta}(1)} = {{C_{S,\theta}(1)} = {- \frac{1}{4}}}},{{C_{C}(1)} = {{S_{S}(1)} = \frac{1}{2}}}} & {{Eq}.\mspace{14mu} 11} \end{matrix}$

Therefore, using the cosine and sine components of the vertical and horizontal forces corresponding to the first-order term of the rotational angle, the cosine component U_(cos)=mr·cos(θ₀) and the sine component U_(sin)=mr·sin(θ₀) of the unbalanced amount can be expressed by the following equation:

$\begin{matrix} {\begin{bmatrix} {C_{F_{x}}(1)} \\ {S_{F_{x}}(1)} \\ {C_{F_{y}}(1)} \\ {S_{F_{y}}(1)} \end{bmatrix} = {\quad{\begin{bmatrix} {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} & {{b \cdot {C_{C}(1)}} - {2\;{b \cdot {S_{C,\theta}(1)}}}} \\ {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} & {{{- a^{2}} \cdot {S_{S}(1)}} - {2\;{b \cdot {S_{S,\theta}(1)}}}} \\ {{{- b} \cdot {C_{C}(1)}} + {2\;{b \cdot {S_{C,\theta}(1)}}}} & {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} \\ {{a^{2} \cdot {S_{S}(1)}} + {2\;{b \cdot {S_{S,\theta}(1)}}}} & {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}}} & {{Eq}.\mspace{14mu} 12} \end{matrix}$

Since an initial angular velocity a and an angular acceleration b are estimated upon every revolution, the angular position θ₀ and the unbalanced amount mr are estimated from the four kinds of force components measured from the support. If two kinds of horizontal components or two kinds of vertical components, or one horizontal component or one vertical component are measured, the balancing is realized. This is obvious that the two unknown quantities, i.e., the cosine and sine components of the unbalanced amount are calculated by only two kinds of measured force components. The conventional balancing method which is performed under the condition of a constant rotating speed is corresponding to a specific case that the angular acceleration is zero, i.e., b=0. Therefore, it can be understood that the balancing method using a linear angular velocity, i.e., a linear time-varying angular velocity model is a further improved method comparing to the conventional method.

In case of higher orders than the number of revolutions, i.e., k≧2, the cosine and sine components of the horizontal and vertical force can be expressed by the follow equations.

$\begin{matrix} {\begin{matrix} {{C_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{{mr} \cdot 2}\;{b \cdot {C_{C,\theta}(k)} \cdot \cos}\left( \theta_{0} \right)} -}} \\ {{{mr} \cdot 2}\;{b \cdot {S_{C,\theta}(k)} \cdot {\sin\left( \theta_{0} \right)}}} \\ {= {{{- {mr}} \cdot 2}\;{b \cdot {S_{C,\theta}(k)} \cdot {\sin\left( \theta_{0} \right)}}}} \end{matrix}\begin{matrix} {{S_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{{+ {mr}} \cdot 2}\;{b \cdot {C_{S,\theta}(k)} \cdot \cos}\left( \theta_{0} \right)} -}} \\ {{{mr} \cdot 2}\;{b \cdot {S_{S,\theta}(k)} \cdot {\sin\left( \theta_{0} \right)}}} \\ {= {{{mr} \cdot 2}\;{b \cdot {C_{S,\theta}(k)} \cdot {\cos\left( \theta_{0} \right)}}}} \end{matrix}\begin{matrix} {{C_{F_{Y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{{mr} \cdot 2}\;{b \cdot {S_{C,\theta}(k)} \cdot {\cos\left( \theta_{0} \right)}}} +}} \\ {{{mr} \cdot 2}\;{b \cdot {C_{C,\theta}(k)} \cdot {\sin\left( \theta_{0} \right)}}} \\ {= {{{mr} \cdot 2}\;{b \cdot {S_{C,\theta}(k)} \cdot {\cos\left( \theta_{0} \right)}}}} \end{matrix}\begin{matrix} {{S_{F_{Y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}} \\ {= {{{{mr} \cdot 2}\;{b \cdot {S_{S,\theta}(k)} \cdot {\cos\left( \theta_{0} \right)}}} +}} \\ {{{mr} \cdot 2}\;{b \cdot {C_{S,\theta}(k)} \cdot {\sin\left( \theta_{0} \right)}}} \\ {= {{{mr} \cdot 2}\;{b \cdot {C_{S,\theta}(k)} \cdot {\sin\left( \theta_{0} \right)}}}} \end{matrix}} & {{Eq}.\mspace{14mu} 13} \end{matrix}$

Like in the equation 5, the values related to the Fourier coefficients are as follows.

$\begin{matrix} {{{{C_{C,\theta}(k)} = {{S_{S,\theta}(k)} = 0}},{{{for}\mspace{14mu} k} = 2},3,\ldots}{{{S_{C,\theta}(k)} = \frac{1}{k^{2} - 1}},{{C_{S,\theta}(k)} = \frac{- k}{k^{2} - 1}}}} & {{Eq}.\mspace{14mu} 14} \end{matrix}$

Therefore, using the cosine and sine components of the vertical and horizontal forces corresponding to higher order terms than the first one of the rotational angle, the cosine component U_(cos)=mr·cos(θ₀) and the sine component U_(sin)=mr·sin(θ₀) of the unbalanced amount can be expressed by the following equation.

$\begin{matrix} {{\begin{bmatrix} {C_{F_{x}}(k)} \\ {S_{F_{x}}(k)} \\ {C_{F_{y}}(k)} \\ {S_{F_{y}}(k)} \end{bmatrix} = {\begin{bmatrix} 0 & {{- 2}\;{b \cdot {S_{C,\theta}(0)}}} \\ {2\;{b \cdot {C_{S,\theta}(k)}}} & 0 \\ {2\;{b \cdot {S_{C,\theta}(2)}}} & 0 \\ 0 & {2\;{b \cdot {C_{S,\theta}(k)}}} \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}},{{{for}\mspace{14mu} k} = 1},2,{3\ldots}} & {{Eq}.\mspace{14mu} 15} \end{matrix}$

Since the angular acceleration b is estimated upon every revolution, the position θ₀ and the unbalanced amount mr are estimated from only two components of the four kinds of force components measured at the support. In other words, if two kinds of horizontal components or two kinds of vertical components, or one horizontal component or one vertical component are measured, the balancing is realized. However, since the conventional balancing method which is performed under the condition of a constant rotating speed is corresponding to the case that the angular acceleration is zero, i.e., b=0, it is not possible to realize the balancing using the cosine and sine components of the horizontal and vertical forces which is corresponding to a second or more-order terms. Therefore, it can be understood that the balancing method using a linear angular velocity, i.e., a linear time-varying angular velocity model is a further improved method comparing to the conventional method.

Optimal Method of Estimating Angle and Magnitude of Unbalanced Mass

Until now, the balancing method using only each separated-order of the zeroth, first and second or more-orders with respect to the rotational angle is described.

Other method of estimating the position θ₀ and the unbalanced amount mr using various-orders with respect to the rotational angle, i.e., k=0, 1, 2, . . . , is described below.

The cosine and sine components of the horizontal and vertical forces measured at the support, and the cosine component U_(cos)=mr·cos(θ₀) and the sine component U_(sin)=mr·sin(θ₀) of the unbalanced amount mr can be expressed by the following equation.

$\begin{matrix} {\begin{bmatrix} {C_{F_{X}}(0)} \\ {C_{F_{Y}}(0)} \\ {C_{F_{X}}(1)} \\ {S_{F_{X}}(1)} \\ {C_{F_{Y}}(1)} \\ {S_{F_{Y}}(1)} \\ {C_{F_{X}}(2)} \\ {S_{F_{X}}(2)} \\ {C_{F_{Y}}(2)} \\ {S_{F_{Y}}(1)} \\ \bullet \\ \bullet \end{bmatrix} = {\quad{\begin{bmatrix} 0 & {{- 2}\;{b \cdot {S_{C,\theta}(0)}}} \\ {{+ 2}\;{b \cdot {S_{C,\theta}(0)}}} & 0 \\ {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} & {{b \cdot {C_{C}(1)}} - {2\;{b \cdot {S_{C,\theta}(1)}}}} \\ {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} & {{{- a^{2}} \cdot {S_{S}(1)}} - {2\;{b \cdot {S_{S,\theta}(1)}}}} \\ {{{- b} \cdot {C_{C}(1)}} + {2\;{b \cdot {S_{C,\theta}(1)}}}} & {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} \\ {{a^{2} \cdot {S_{S}(1)}} + {2\;{b \cdot {S_{S,\theta}(1)}}}} & {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} \\ 0 & {{- 2}\;{b \cdot {S_{C,\theta}(2)}}} \\ {2\;{b \cdot {C_{S,\theta}(2)}}} & 0 \\ {2\;{b \cdot {S_{S,\theta}(2)}}} & 0 \\ 0 & {2\;{b \cdot {C_{S,\theta}{()}}}} \\ \bullet & \bullet \\ \bullet & \bullet \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}}} & {{Eq}.\mspace{14mu} 16} \end{matrix}$

The left member of the equation indicates cosine and sine items by orders of the Fourier coefficients of the horizontal and vertical forces applied to the support. An unbalanced vector of the right member, which has the size of 2-by-1, indicates the cosine and sine components of the unbalanced amount and can be expressed as follows.

$\begin{matrix} {{U_{\cos} = {{U \cdot {\cos\left( \theta_{0} \right)}} = {{mr} \cdot {\cos\left( \theta_{0} \right)}}}},{U_{\sin} = {{U \cdot {\sin\left( \theta_{0} \right)}} = {{mr} \cdot {\sin\left( \theta_{0} \right)}}}}} & {{Eq}.\mspace{14mu} 17} \end{matrix}$

In the equation 10, assuming that the magnitude of left vector V_(F) is L-by-1, the right matrix M has the size of L-by-2 and can be expressed by the following matrix.

$\begin{matrix} {V_{F} = {M \cdot \begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}} & {{Eq}.\mspace{14mu} 18} \end{matrix}$

The matrix of the equation 18 represents a linear relationship between each order component of the Fourier coefficients of the horizontal and vertical forces measured at the support 30 and the two unknown quantities, i.e., the cosine and sine components of the unbalanced amount. From a plurality of measured values having the linear relationship, the two unknown optimal quantities are determined by the least squares method. The following equation which is expressed by a determinant is to calculate the two unknown optimal quantities.

$\begin{matrix} {\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix} = {M^{+} \cdot V_{F}}} & {{Eq}.\mspace{14mu} 19} \end{matrix}$

In the equation 19, M⁺ indicates a pseudo inverse matrix for the matrix M. The calculation of the pseudo inverse matrix is readily performed according to the numerical method well-known in linear algebra. The unbalanced amount and the angular position are optimally estimated from the cosine and sine components of the unbalanced mass obtained by the series of calculations, as follows.

$\begin{matrix} {{U = {{mr} = \sqrt{U_{\cos}^{2} + U_{\sin}^{2}}}},{\theta_{0} = {\tan^{- 1}\left( \frac{U_{\sin}}{U_{\cos}} \right)}}} & {{Eq}.\mspace{14mu} 20} \end{matrix}$

The above result is an optimal method of estimating the unbalanced amount and the angular position by measuring the components of the horizontal and vertical forces applied to the support. Most of the balancing apparatuses used on the ground estimate the unbalanced amount and the angular position by measuring only one component corresponding to the horizontal and vertical forces. To estimate the unbalanced amount and the angular position using only one of the horizontal and vertical forces, the values corresponding to the non-measured terms in the general equation are substituted by zero, or the terms are eliminated, and then the least square method or the above-mentioned determinant is performed or calculated. The following equation indicates the case that only the horizontal force component is measured.

$\begin{matrix} {\begin{bmatrix} {C_{F_{X}}(0)} \\ {C_{F_{X}}(1)} \\ {S_{F_{X}}(1)} \\ {C_{F_{X}}(2)} \\ {S_{F_{X}}(2)} \\ \bullet \end{bmatrix} = {\quad{\begin{bmatrix} 0 & {{- 2}\;{b \cdot {S_{C,\theta}(0)}}} \\ {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} & {{b \cdot {C_{C}(1)}} - {2\;{b \cdot {C_{C,\theta}(1)}}}} \\ {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} & {{{- a^{2}} \cdot {S_{S}(1)}} - {2\;{{bC}_{S,\theta}(1)}}} \\ 0 & {{- 2}\;{b \cdot {S_{C,\theta}(0)}}} \\ {2\;{b \cdot {C_{S,\theta}(2)}}} & 0 \\ \bullet & \bullet \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}}} & {{Eq}.\mspace{14mu} 21} \end{matrix}$

As shown in the equation 21, the cosine and sine terms have the linear relationship with the various order terms of the Fourier coefficients of the horizontal force. Therefore, it can be understood that the optimal cosine and sine components of the unbalanced mass are estimated from a single directional force. Also, the angle θ₀ and the magnitude U of the unbalanced amount are calculated by the above-mentioned equation. In case that only the vertical force is measured, it is possible to estimate the optimal co sine and sine components of the unbalanced mass by the same method.

Method of Sampling Applied Force Signal and Method of Estimating Angular Velocity Coefficient

If a rotary encoder is installed at a rotational shaft of the balancing apparatus, it is possible to obtain a constant angle signal upon every revolution of the shaft. In other words, it is possible to obtain N equi-angular pulse signals per revolution of the shaft, i.e., an angular position pulse having a constant width along a rotational angle. The angular position pulses are used as a synchronous pulse for sampling the applied force signals of the support at every desired angle. The method of sampling the signal at each of N equi-angles per revolution is called the equi-angle sampling method.

Hereinafter, there is introduced a method of estimating the angular velocity and the angular acceleration at every revolution of the rotational shaft using a reference position pulse generated whenever the rotary encoder is rotated once. Assuming that each revolution period before and after the reference position pulse which is indicated by “0” in FIG. 2 is respectively T_(n−1)(T_(n−1)=N_(n−1)·ΔT) and T_(n)(T_(n)=N_(n)·ΔT), an initial angular velocity a_(n) and an angular acceleration b_(n) can be expressed by the following equation.

$\begin{matrix} {{a_{n} = {\frac{2\;\pi}{\Delta\; T} \cdot \frac{N_{n - 1}^{2} + N_{n}^{2}}{N_{n - 1}{N_{n} \cdot \left( {N_{n - 1} + N_{n}} \right)}}}},{b_{n} = {\frac{4\;\pi}{\Delta\; T^{2}} \cdot \frac{N_{n - 1} - N_{n}}{N_{n - 1}{N_{n} \cdot \left( {N_{n - 1} + N_{n}} \right)}}}}} & {{Eq}.\mspace{14mu} 22} \end{matrix}$

In the equation 22, ΔT indicates a period of a base clock used to measure a period of the reference position pulse, and mainly uses a frequency of 20 MHz or more. The initial angular velocity a_(n) and the angular acceleration b_(n) are respectively corresponding to the initial angular velocity a and the angular acceleration b in the linear time-varying angular velocity model.

If the rotary encoder is not used, the applied force signal of the support is sampled by a constant time sampling method which is employed in most of digital signal sampling methods. As shown in FIG. 3, the constant time sampling method uses a sampling pulse having a constant time width (Δt=constant). In case of sampling the applied force output signal to be measured with a constant time width, as shown in FIG. 3, it is substantially impossible to synchronize the reference position pulse and the sampling pulse. Therefore, the present invention proposes a method of calculating a measurement value at a reference point by measuring a time interval between the reference position pulse and the sampling pulse, and also calculating an angle between the reference point and a first sampling pulse.

In the lower right side of FIG. 3, a time interval d_(n)·Δt between an n-th reference position pulse and a sampling pulse is illustrated, and in the lower left side of FIG. 3, an equation model for calculating a measured force component F_(X,0) corresponding to a time when the n-th reference position pulse is generated, using the interpolation principle, is provided. A time interval between the reference position pulse and the first sampling signal of every revolution is measured by a period measuring circuit of the reference position pulse as shown in FIG. 4.

In FIG. 4, C_(n) and D_(n) indicate a period of the reference pulse and a digital value used in measuring an observing time of a first sampling pulse after the reference pulse, respectively. The period of the reference position pulse and the time difference between the reference pulse and the first sampling pulse are calculated by the following equation: T _(n)=(C _(n+1) −C _(n))·ΔT=N _(n) ·ΔT, d _(n) ·Δt=(D _(n) −C _(n))·ΔT  Eq. 23

Using the initial angular velocity a_(n) and the angular acceleration b_(n) corresponding to the n-th revolution period and the time interval d_(n)·Δt between the n-th reference position pulse and the sampling signal, the rotational angle between the sampling pulses can be expressed as follows: Δθ₀=θ(d _(n) ·Δt), Δθ(m·Δt)=θ((m+1)·Δt)−θ(m·Δt), m=1, 2, . . . , M _(n)−1 Δθ_(M) _(n) =2π−θ((d _(n) +M _(n)−1)·Δt).  Eq. 24

In the equation 24, M_(n) is the number of signals sampled during an n-th revolution period. By using the rotational difference between the sampling pulse, the cosine value of the Fourier coefficient of the horizontal force components measured at the support is calculated.

$\begin{matrix} {{C_{F_{X}}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \ {F_{X}(\theta)}}{\mathbb{d}\theta}}}} \cong {\frac{1}{2\;\pi} \cdot \left\{ {{{F_{X,0} \cdot \Delta}\;\theta_{0}} + {\sum\limits_{m = 1}^{M_{n} - 1}{{{F_{X}\left( {{m \cdot \Delta}\; t} \right)} \cdot \Delta}\;\theta\left( {{m \cdot \Delta}\; t} \right)}} + {{{F_{X}\left( {{M_{n} \cdot \Delta}\; t} \right)} \cdot \Delta}\;\theta_{0}}} \right\}}}} & {{Eq}.\mspace{14mu} 25} \end{matrix}$

In the equation 25, F_(X,0) is the horizontal force component corresponding to a time when the reference position pulse passes, and is calculated by using the interpolation principle as illustrated in the lower left side of FIG. 3. Of course, the cosine and sine values of the Fourier coefficients of the vertical force components as well as the sine value of the Fourier coefficients of the horizontal force components are estimated by the numerical integration method which is the same as the method of calculating the cosine value of the Fourier coefficients of the horizontal force components shown above. To reduce an error between a numerical integration value and a theoretical value, it is preferable that a sampling rate is established such that the number M_(n) of the signals sampled during one revolution is at least 100 or more.

Balancing Using Vibration Measurement of Support

Until now, the method of directly measuring the applied forces on the support 30 to have the condition of a rigid body is described. If the supporting condition of the rotational shaft is minutely deformed or vibrated like a hard-spring having sufficient rigidity, the linear relationship is maintained between the force applied to the support and the vibration of the support. FIGS. 5A, 5B and 5C show three kinds of vibration models of the support, i.e., a horizontal vibration model, a vertical vibration model and a dual-axis vibration model.

The single direction models, i.e., the horizontal vibration model and the vertical vibration model are most widely used in the horizontal and vertical balancing apparatus, and the dual-axis vibration model is also used in the precision vertical balancing apparatus, i.e., a model which needs horizontal direction rigidity. As known in the field of vibration, the vibration acceleration and the force applied to the support have a linear proportion relationship, as follows. F _(X) ≈m _(eff,X) ·a _(X)(t), F _(Y) ≈m _(eff,Y) ·a _(Y)(t)  Eq. 26

In the equation 26, m_(eff) indicates an effective mass of each measured shaft, and the effective mass is regarded as a proportional constant representing the linear relationship between the magnitude of the acceleration measured and the applied force to the practical test piece. This process is called the field balancing calibration using a known test mass and when the calibration is completed, the vibration acceleration component measured by a vibration sensor is readily used to evaluate the force applied to the support. Therefore, the methods of estimating the angular position and the magnitude of the unbalanced amount using the above-mentioned measured force are applied to the estimation of the angular position and the unbalanced amount of the rotational body using the Fourier coefficients of an output signal of the vibration sensor after the completion of the field balancing calibration. In other words, the method of estimating the angle and unbalanced amount, in which periodic rotational vibration is applied to the test piece and the single or dual axis vibration sensors are installed at the left and right (or upper and lower) support having the desired rigidity, is carried out along the same procedure as in the estimating method by directly measuring the applied forces.

Balancing Apparatus

FIG. 6 is a schematic view showing the configuration of a dynamic balancing apparatus using a linear time-varying angular velocity model in accordance with an embodiment of the present invention.

As described in FIG. 6, the dynamic balancing apparatus using a linear time-varying angular velocity model of the present invention includes a rotational shaft 10 on which a rotational body 1 having an unbalanced mass is installed; a linear time-varying angular velocity generator 20 which allow the rotational shaft 10 to be rotated at a linear time-varying angular velocity; a support 30 for supporting both sides of the rotational shaft 10; and a transducer 40 for measuring applied force or vibration transferred to the support 30 by revolution of the rotational shaft 10 having a linear time-varying angular velocity.

The rotational shaft 10 has the rotational body 1 having an unbalanced mass. The linear time-varying angular velocity generator 20 is installed at the rotational shaft 10 so that the rotational body 1 on the rotational shaft 10 is rotated at the linear time-varying angular velocity.

The linear time-varying angular velocity generator 20 functions to generate the rotational vibration due to the rotating unbalanced mass of the rotational body 1 so that the rotational shaft 10 is rotated at the linear time-varying angular velocity.

The support 30 supports the both sides of the rotational shaft 10. The rotational vibration generated by the linear time-varying angular velocity or the applied forces, i.e., centrifugal and tangential forces, are transferred to the support 30 by the revolution of the rotational shaft 10.

The transducer 40 is installed at the support 30 so as to measure the applied force or vibration transferred to the support 30 by the revolution or vibration of the rotational shaft 10. Thus, it is possible to estimate the angle and magnitude of the rotating unbalanced mass of the rotational body 1 using characteristic of the applied force or vibration measured by the transducer 40.

In addition, it is preferable that the present invention further includes a reference pulse generator 50 which generates a reference position pulse of the rotational shaft 10 so as to measure a position of the rotational body 1 with respect to the linear time-varying angular velocity of the rotational shaft 10.

Unlike the conventional balancing apparatus, the dynamic balancing apparatus using a linear time-varying angular velocity model uses the linear time-varying angular velocity generator 20 and a rotational vibration measuring sensor for measuring the amplitude of the periodic rotational vibration.

Angular motion controllers mainly allow a user to precisely control the linear time-varying angular velocity. Since high-speed motors like a brushless direct current (BLDC) motor or high-torque motors mainly provide a wide torque range of a few tens rpm to a few thousands rpm, these motors may be used as an apparatus for generating the linear angular motion. In case of using the reference pulse generator 50 or the rotary encoder which measures the reference position of the rotational shaft and the position of the rotational body at each constant angle, it is possible to sample the certain number of signals at every revolution with respect to the components of the horizontal or vertical force applied to the support using the equi-angle sampling method as described above. In case that the rotary encoder or the reference pulse generator 50 is not installed at the balancing apparatus, an apparatus for generating one reference position pulse upon every revolution is provided at the rotational shaft, and then the balancing is carried out using the constant time sampling method.

The apparatus for generating the linear angular velocity motion, the rotary encoder which is installed at the rotational shaft or, the apparatus for generating one reference position pulse may be additionally provided at the existing balancing apparatus so as to carry out the balancing. This means that the existing balancing apparatus can also perform the balancing method proposed in the present invention by appropriately applying the method of controlling the linear angular velocity motion with respect to time.

According to the present invention, it is possible to solve the problem that the conventional balancing apparatus performs the balancing only in the state that a motor is continuously rotated at a constant velocity. Therefore, since the present invention can perform the balancing even in the accelerating time until reaching a certain rotating speed from the stopped state and the decelerating time which is necessary to stop the rotating state from the certain rotating speed, it can be applied to an automatic balancing apparatus for a rapid production process. Thus, it is possible to remarkably reduce a balancing time in the automatic balancing apparatus for rapid production processes.

Those skilled in the art will appreciate that the conceptions and specific embodiments disclosed in the foregoing description may be readily utilized as a basis for modifying or designing other embodiments for carrying out the same purposes of the present invention. Those skilled in the art will also appreciate that such equivalent embodiments do not depart from the spirit and scope of the invention of the appended claims. 

The invention claimed is:
 1. A dynamic balancing method using a linear time-varying angular velocity model, comprising: rotating a rotational shaft, on which a rotational body having a rotating unbalanced mass is installed, at a linear time-varying angular velocity; measuring vibration generated by tangential force, centrifugal force, or the two forces generated by the rotating unbalanced mass of the rotational body from the support which supports the rotational shaft rotated at the linear time-varying angular velocity; estimating multiple frequency components of a measurement signal corresponding to the measured normal and centrifugal forces; and estimating the angular position, that is the angle, and the magnitude of the unbalanced mass of the rotational body using components of the measured tangential and centrifugal forces and characteristics of the frequency components (or Fourier coefficients) of the two forces, wherein the angular position and magnitude of the unbalanced mass of the rotational body is determined by the following equation: $\begin{matrix} {{F_{X}\left( {\theta(t)} \right)} = {{{- F_{t}} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} + {F_{n} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}}}} \\ {= {{{mr} \cdot b \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} + {{mr} \cdot \left\{ {a^{2} + {2\;{b \cdot {\theta(t)}}}} \right\} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}}}} \\ {= {{mr} \cdot b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\cos\left( {\theta(t)} \right\}}} +} \right.}} \right.}} \\ {{mr} \cdot a^{2} \cdot \left\{ {{{\cos\left( \theta_{0} \right)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\sin(\theta)}}} \right\}}} +} \right.} \\ {{{mr} \cdot 2}\;{b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\sin\left( {\theta(t)} \right\}}}} \right.}} \right.}} \end{matrix}$ $\begin{matrix} {{F_{Y}\left( {\theta(t)} \right)} = {{F_{t} \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}} + {F_{n} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} -}} \\ {= {{{mr} \cdot b \cdot {\cos\left( {{\theta(t)} + \theta_{0}} \right)}} + {{mr} \cdot \left\{ {a^{2} + {2\;{b \cdot {\theta(t)}}}} \right\} \cdot {\sin\left( {{\theta(t)} + \theta_{0}} \right)}} -}} \\ {= {{mr} \cdot b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\cos\left( {{\theta(t)} - {{\sin\left( \theta_{0} \right)} \cdot {\sin\left( {\theta(t)} \right\}}} +} \right.}} \right.}} \\ {{mr} \cdot a^{2} \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\cos\left( {\theta(t)} \right\}}} +} \right.}} \right.} \\ {{{mr} \cdot 2}\;{b \cdot \left\{ {{\cos\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\sin\left( {{\theta(t)} + {{\sin\left( \theta_{0} \right)} \cdot {\theta(t)} \cdot {\cos\left( {\theta(t)} \right\}}}} \right.}} \right.}} \end{matrix}$ wherein F_(x)(t) is the horizontal component of each of the tangential force F_(t) and the centrifugal force F_(n) acting on the support by the rotating unbalanced mass of the rotational body; F_(y)(t) is the vertical component of each of the tangential force F_(t) and the centrifugal force F_(n) generated acting on the support by the rotating unbalanced mass of the rotational body; m[kg] is the rotating unbalanced mass of the rotational body; r[m] is a distance of the rotating unbalanced mass of the rotational body, which is spaced apart from the center of rotational shaft; a[radians/s] is an initial rotational velocity; b[radians/s²] is a constant angular acceleration at the initial rotational velocity; θ(t) is an angular displacement of the rotational shaft; and θ₀[radians] is an angular position from the rotation reference point of the rotational shaft wherein the angular position and magnitude of the unbalanced mass of the rotational body is determined by the following equation using the tangential and centrifugal forces, which are measured in the support and applied in the horizontal and vertical directions, ${\begin{bmatrix} {C_{F_{x}}(k)} \\ {S_{F_{x}}(k)} \\ {C_{F_{y}}(k)} \\ {S_{F_{y}}(k)} \end{bmatrix} = {\begin{bmatrix} 0 & {{- 2}\;{b \cdot {S_{C,\theta}(0)}}} \\ {2\;{b \cdot {C_{S,\theta}(k)}}} & 0 \\ {2\;{b \cdot {S_{C,\theta}(2)}}} & 0 \\ 0 & {2\;{b \cdot {C_{S,\theta}(k)}}} \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}},{{{for}\mspace{14mu} k} = 1},2,{3\ldots}$ wherein C_(Fx) is the cosine component corresponding to Fourier coefficients of the horizontal component of the tangential force and the centrifugal force acting on the support by the rotating unbalanced mass of the rotational body, and S_(Fx) is the sine component corresponding to Fourier coefficients of the horizontal component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and these components are determined by the following equation: ${{C_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{X}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and k is an integer indicating the order of Fourier coefficients, C_(Fy) is the cosine component corresponding to Fourier coefficients of the vertical component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and S_(Fy) is the sine component corresponding to Fourier coefficients of vertical components of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and these components are determined by the following equation: ${{C_{F_{y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and U_(cos) is the cosine component of the unbalanced amount (U_(cos)=mr·cos(θ₀)) and U_(sin) is the sine component of the unbalanced amount (U_(sin)=mr·sin(θ₀)), ${{C_{c,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${{S_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${C_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ {{\begin{matrix} {{{- 1}/4},\mspace{14mu}{{{for}\mspace{14mu} k} = 1}} \\ {{- \frac{k}{k^{2} - 1}},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,\ldots} \end{matrix}{S_{c,\theta}(k)}} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ \begin{matrix} {{{3\;{k/4}} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 0},1} \\ {\frac{1}{k^{2} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,{\ldots.}} \end{matrix} \right.}} \right.}$ wherein the angular position and magnitude of the unbalanced mass of the rotational body is determined by the following equation using the tangential and centrifugal forces applied in the horizontal direction, $\begin{bmatrix} {C_{F_{x}}(0)} \\ {C_{F_{x}}(1)} \\ {S_{F_{x}}(1)} \\ {C_{F_{x}}(2)} \\ {S_{F_{x}}(2)} \\ \bullet \end{bmatrix} = {\begin{bmatrix} 0 & {{- 2}{b \cdot {S_{C,\theta}(0)}}} \\ {{a^{2} \cdot {C_{c}(1)}} + {2{b \cdot {C_{C,\theta}(1)}}}} & {{b \cdot {C_{C}(1)}} - {2{b \cdot {C_{C,\theta}(1)}}}} \\ {{b \cdot {S_{S}(1)}} + {2{b \cdot {C_{S,\theta}(1)}}}} & {{{- a^{2}} \cdot {S_{S}(1)}} - {2{{bC}_{S,\theta}(1)}}} \\ 0 & {{- 2}{b \cdot {S_{C,\theta}(2)}}} \\ {2{b \cdot {C_{S,\theta}(2)}}} & 0 \\ \bullet & \bullet \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}$ wherein C_(Fx) is the cosine component corresponding to Fourier coefficients of the horizontal component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and S_(Fx) is the sine component corresponding to Fourier coefficients of the horizontal component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and these components are determined by the following equation: ${{C_{F_{x}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{x}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and k is an integer indicating the order of the Fourier coefficient, U_(cos) is the cosine component of the unbalanced amount (U_(cos)=mr·cos(θ₀)) and U_(sin) is the sine component of the unbalanced amount (U_(sin)=mr·sin(θ₀)), and ${{C_{c,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${{S_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${C_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ {{\begin{matrix} {{{- 1}/4},\mspace{14mu}{{{for}\mspace{14mu} k} = 1}} \\ {{- \frac{k}{k^{2} - 1}},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,\ldots} \end{matrix}{S_{c,\theta}(k)}} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ \begin{matrix} {{{3\;{k/4}} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 0},1} \\ {\frac{1}{k^{2} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,{\ldots.}} \end{matrix} \right.}} \right.}$ wherein a measurement signal sampling method exploits a sampling pulse having a constant time width (Δt=constant) and then measures applied force signals of the support using the following equation: T _(n)=(C _(n+1) −C _(n))·ΔT=N _(n) ·ΔT, d _(n) ·Δt=(D _(n) −C _(n))·ΔT wherein C_(n) and D_(n) are a period of the reference pulse and a digital value used in measuring an observing time of a first sampling pulse after the reference pulse, ΔT is a period of a base clock used in measuring a period of the reference position pulse, T_(n) is n-th revolution period, and d_(n) is a fractional part of a sampling period, d_(n) is a fractional part of a sampling period, and d_(n)·Δt is a time interval between an n-th reference position pulse and a sampling pulse, and Δθ₀=θ(d _(n) ·Δt), Δθ(m·Δt)=θ((m+1)·Δt)−θ(m·Δt), m=1, 2, . . . , M _(n)−1 Δθ_(M) _(n) =2π−θ((d _(n) +M _(n)−1)·Δt) M_(n) is the number of signals sampled during one revolution, and ${C_{F_{x}}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \ {F_{X}(\theta)}}{\mathbb{d}\theta}}}} \cong {\frac{1}{2\;\pi} \cdot \begin{Bmatrix} {{{F_{X,0} \cdot \Delta}\;\theta_{0}} + {\sum\limits_{m = 1}^{M_{n} - 1}{{F_{X}\left( {{m \cdot \Delta}\; t} \right)} \cdot}}} \\ {{\Delta\;\theta\left( {{m \cdot \Delta}\; t} \right)} + {{{F_{X}\left( {{M_{n} \cdot \Delta}\; t} \right)} \cdot \Delta}\;\theta_{0}}} \end{Bmatrix}}}$ C_(Fx)(k) is the cosine value of Fourier coefficients of the component of the horizontal force measured in the support.
 2. The dynamic balancing method of claim 1, wherein the angular position and magnitude of the unbalanced mass of the rotational body is determined by the following equation using the tangential force and the centrifugal force acting on support in the horizontal and vertical directions, $\begin{bmatrix} {C_{F_{X}}(0)} \\ {C_{F_{Y}}(0)} \\ {C_{F_{X}}(1)} \\ {S_{F_{X}}(1)} \\ {C_{F_{Y}}(1)} \\ {S_{F_{Y}}(1)} \\ {C_{F_{X}}(2)} \\ {S_{F_{X}}(2)} \\ {C_{F_{Y}}(2)} \\ {S_{F_{Y}}(1)} \\ \begin{matrix} \bullet \\ \bullet \end{matrix} \end{bmatrix} = {\begin{bmatrix} 0 & {{- 2}\;{b \cdot {S_{C,\theta}(0)}}} \\ {{+ 2}\;{b \cdot {S_{C,\theta}(0)}}} & 0 \\ {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} & {{b \cdot {C_{C}(1)}} - {2\;{b \cdot {S_{C,\theta}(1)}}}} \\ {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} & {{{- a^{2}} \cdot {S_{S}(1)}} - {2\;{b \cdot {S_{S,\theta}(1)}}}} \\ {{{- b} \cdot {C_{C}(1)}} + {2\;{b \cdot {S_{C,\theta}(1)}}}} & {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} \\ {{a^{2} \cdot {S_{S}(1)}} + {2\;{b \cdot {S_{S,\theta}(1)}}}} & {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} \\ 0 & {{- 2}\;{b \cdot {S_{C,\theta}(2)}}} \\ {2\;{b \cdot {C_{S,\theta}(2)}}} & 0 \\ {2\;{b \cdot {S_{S,\theta}(2)}}} & 0 \\ 0 & {2\;{b \cdot {C_{S,\theta}{()}}}} \\ \bullet & \bullet \\ \bullet & \bullet \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}$ wherein C_(FX) is the cosine component corresponding to Fourier coefficients of the horizontal component of the tangential force and the centrifugal force acting on the support by the rotating unbalanced mass of the rotational body, and S_(Fx) is the sine component corresponding to Fourier coefficients of the horizontal component of the tangential and centrifugal forces generated at the support by the rotating unbalanced mass of the rotational body, and these components are determined by the following equation: ${{C_{F_{x}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{x}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{X}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and k is an integer indicating the order of Fourier coefficients, C_(Fy) is the cosine component corresponding to Fourier coefficients of the vertical component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and S_(Fy) is the sine component corresponding to Fourier coefficients of the vertical component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and these components are determined by the following equation: ${{C_{F_{y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and U_(cos) is the cosine component of the unbalanced amount (U_(cos)=mr·cos(θ₀)) and U_(sin) is the sine component of the unbalanced amount (U_(sin)=mr·sin(θ₀)), and ${{C_{c,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${{S_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${C_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ {{\begin{matrix} {{{- 1}/4},\mspace{14mu}{{{for}\mspace{14mu} k} = 1}} \\ {{- \frac{k}{k^{2} - 1}},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,\ldots} \end{matrix}{S_{c,\theta}(k)}} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ \begin{matrix} {{{3\;{k/4}} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 0},1} \\ {\frac{1}{k^{2} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,{\ldots.}} \end{matrix} \right.}} \right.}$
 3. The dynamic balancing method of claim 1, wherein the angular position and magnitude of the unbalanced mass of the rotational body is determined by the following equation using the tangential force and the centrifugal force applied in the vertical direction, $\begin{bmatrix} {C_{F_{y}}(0)} \\ {C_{F_{y}}(1)} \\ {S_{F_{y}}(1)} \\ {C_{F_{y}}(2)} \\ {S_{F_{y}}(1)} \\ \bullet \\ \bullet \end{bmatrix} = {\begin{bmatrix} {{+ 2}\;{b \cdot {S_{C,\theta}(0)}}} & 0 \\ {{{- b} \cdot {C_{C}(1)}} + {2\;{b \cdot {S_{C,\theta}(1)}}}} & {{a^{2} \cdot {C_{C}(1)}} + {2\;{b \cdot {C_{C,\theta}(1)}}}} \\ {{a^{2} \cdot {S_{S}(1)}} + {2\;{b \cdot {S_{S,\theta}(1)}}}} & {{b \cdot {S_{S}(1)}} + {2\;{b \cdot {C_{S,\theta}(1)}}}} \\ {2\;{b \cdot {S_{S,\theta}(2)}}} & 0 \\ 0 & {2\;{b \cdot {C_{S,\theta}(2)}}} \\ \bullet & \bullet \\ \bullet & \bullet \end{bmatrix}\begin{bmatrix} U_{\cos} \\ U_{\sin} \end{bmatrix}}$ wherein C_(Fy) is the cosine component corresponding to Fourier coefficients of the vertical component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and S_(Fy) is the sine component corresponding to Fourier coefficients of the vertical component of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body, and these components are determined by the following equation: ${{C_{F_{y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}},{{S_{F_{y}}(k)} = {\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot {F_{Y}(\theta)}}\ {\mathbb{d}\theta}}}}}$ and k is an integer indicating the order of the Fourier coefficient, and U_(cos) is the cosine component of the unbalanced amount (U_(cos)=mr·cos(θ₀)) and U_(sin) is the sine component of the unbalanced amount (U_(sin)=mr·sin(θ₀)), and ${{C_{c,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${{S_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \frac{\pi}{2}}},\mspace{14mu}{{{only}\mspace{14mu}{for}\mspace{14mu} k} = 1}$ ${C_{s,\theta}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\sin\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\cos(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ {{\begin{matrix} {{{- 1}/4},\mspace{14mu}{{{for}\mspace{14mu} k} = 1}} \\ {{- \frac{k}{k^{2} - 1}},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,\ldots} \end{matrix}{S_{c,\theta}(k)}} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \theta \cdot {\sin(\theta)}}\ {\mathbb{d}\theta}}}} = \left\{ \begin{matrix} {{{3\;{k/4}} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 0},1} \\ {\frac{1}{k^{2} - 1},\mspace{14mu}{{{for}\mspace{14mu} k} = 2},3,{\ldots.}} \end{matrix} \right.}} \right.}$
 4. The dynamic balancing method of claim 1, wherein a measurement signal sampling method is used to sample force signals applied to the support using a reference pulse generator or a rotary encoder which generates equi-angle position pulses having a constant angular revolution using a reference position pulse generated whenever the rotational shaft is rotated once.
 5. The dynamic balancing method of claim 4, wherein an angular velocity and an angular acceleration corresponding to a revolution period when generating N equi-angular pulses per revolution are determined by the following equation: ${a_{n} = {\frac{2\;\pi}{\Delta\; T} \cdot \frac{N_{n - 1}^{2} + N_{n}^{2}}{N_{n - 1}{N_{n} \cdot \left( {N_{n - 1} + N_{n}} \right)}}}},{b_{n} = {\frac{4\;\pi}{\Delta\; T^{2}} \cdot \frac{N_{n - 1} - N_{n}}{N_{n - 1}{N_{n} \cdot \left( {N_{n - 1} + N_{n}} \right)}}}}$ wherein a_(n) is an angular velocity corresponding to the revolution period T_(n) of the rotational shaft, b_(n) is an angular acceleration corresponding to the revolution period T_(n) of the rotational shaft, ΔT is a period of a base clock used in measuring a period of the reference position pulse, N_(n−1) is the number of equi-angular pulse signals corresponding to a revolution period T_(n−1)(T_(n−1)=N_(n−1)·ΔT) before the reference position pulse, and N_(n) is the number of equi-angle pulses each revolution period T_(n)(T_(n)=N_(n)·ΔT) after the reference position pulse.
 6. The dynamic balancing method of claim 1, wherein horizontal and vertical components of the tangential and centrifugal forces acting on the support by the rotating unbalanced mass of the rotational body using vibration acceleration measured from a vibration sensor is determined by the following equation: F _(X) ≈m _(eff,X) ·a _(X)(t), F _(Y) ≈m _(eff,Y) ·a _(Y)(t) m_(eff) is an effective mass of a rotational body, and a_(x)(t) and a_(y)(t) are the measured vibration acceleration signals that satisfy a linear relationship between the measured vibration accelerations and the tangential and centrifugal forces acting on the support and the measured vibration acceleration signals from the support are used to determine the angular position and magnitude of an unbalanced mass, by installing vibration transducers in the horizontal and/or vertical movement directions of the support, and Fourier coefficients of the horizontal and vertical force components, are evaluated for implementing vibration measurement-based dynamic balancing.
 7. The dynamic balancing method of claim 3, wherein a measurement signal sampling method is used to sample force signals applied to the support using a reference pulse generator or a rotary encoder which generates equi-angle position pulses having a constant angular revolution using a reference position pulse generated whenever the rotational shaft is rotated once.
 8. The dynamic balancing method of claim 3, wherein a measurement signal sampling method exploits a sampling pulse having a constant time width (Δt=constant) and then measures applied force signals of the support using the following equation: T _(n)=(C _(n+1) −C _(n))·ΔT=N _(n) ·ΔT, d _(n) ·Δt=(D _(n) −C _(n))·ΔT wherein C_(n) and D_(n) are a period of the reference pulse and a digital value used in measuring an observing time of a first sampling pulse after the reference pulse, ΔT is a period of a base clock used in measuring a period of the reference position pulse, T_(n) is n-th revolution period, and d_(n) is a fractional part of a sampling period, and Δθ=θ(d _(n) ·Δt), Δθ=(m·Δt)=θ((m+1)·Δt)−θ(m·Δt), m=1, 2, . . . M _(n)−1 Δθ_(M) _(n) =2π−θ((d _(n) +M _(n)−1)·Δt), M_(n) is the number of signals sampled during one revolution, and ${C_{F_{x}}(k)} = {{\frac{1}{2\;\pi}{\int_{0}^{2\;\pi}{{{\cos\left( {k \cdot \theta} \right)} \cdot \ {F_{X}(\theta)}}{\mathbb{d}\theta}}}} \cong {\frac{1}{2\;\pi} \cdot \begin{Bmatrix} {{{F_{X,0} \cdot \Delta}\;\theta_{0}} + {\sum\limits_{m = 1}^{M_{n} - 1}{{F_{X}\left( {{m \cdot \Delta}\; t} \right)} \cdot}}} \\ {{\Delta\;\theta\left( {{m \cdot \Delta}\; t} \right)} + {{{F_{X}\left( {{M_{n} \cdot \Delta}\; t} \right)} \cdot \Delta}\;\theta_{0}}} \end{Bmatrix}}}$ C_(Fx)(k) is the cosine value of Fourier coefficients of the component of the horizontal force measured in the support. 